Destruction of bulk ordering by surface randomness 
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We demonstrate that the arbitrarily weak quenched disorder on the surface of a system of 
continuous symmetry destroys long range order in the bulk, and, instead, quasi-long range order 
emerges. Correlation functions are calculated exactly for the two- and three-dimensional XY model 
with surface randomness via the functional renormalization group. Even at strong quenched disorder 
the three-dimensional XY model possesses topological order. We also determine roughness of a 
£h ■ domain wall in the presence of surface disorder. 
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The arbitrarily weak quenched disorder in the bulk of a system of continuous symmetry destroys long range order 
inherent to the pure system Q provided that disorder breaks not only the translational symmetry but also the 
symmetry with respect to transformations of the order parameter, as e.g. random anisotropy in amorphous magnets 
does. This fundamental fact governs all the physics of condensed matter and results in a wealth of observed static 

■ and dynamics behaviors of real solids. 

In many cases noticeable disorder presents only at the surface. Not surprisingly, surface randomness modifies the 
critical behavior near the surface Q], yet the common expectation is for the bulk properties to remain intact. In this 
Letter we show that arbitrarily weak surface disorder destroys long range order in the bulk of a system of continuous 

■ symmetry at the arbitrarily low temperature. 

The predicted effect occurs in a rich variety of systems. Examples include crystal ordering in solids grown on a 
*^ . disordered substrate, liquid crystals interacting with an inhomogeneous surface, superconducting vortices pinned by 
surface impurities, etc. There are also many two-dimensional systems with edge randomness, e.g. superconducting 
films with columnar defects in a part of the film or films with a rough edge || . 

The reason as to why surface impurities, however weak, break long range bulk order is that the bulk contributes 
little to the energy of long-wave Goldstone modes: the surface random energy of long-wave excitations turns out to 
be greater than the corresponding bulk energy. As a result, the inhomogeneous state becomes favorable energetically. 
Note that ordering survives in the regions of the size less than the distance of these regions from the surface. In other 
words, if the distance between the two points is greater than their separation from the surface, the order parameter is 
^sO , different in those points. While long range order breaks down, topological order survives and quasi-long range order 
emerges. This means that the correlation length is infinite and that the correlation functions obey a slow logarithmic 
dependence of the distance. 

An easy way to understand the main result of the Letter is based on Imry-Ma arguments Let us consider 

a region of size L near the surface and compare the energies of ordered and disordered states of the region. If long 
range order is broken on the scales of the order of L the loss in the bulk (elastic) energy is -Ebuik ~ L D /L 2 , where D 
is the space dimension. The energy gain from the interaction with surface impurities scales as Surface ~ 
D < 3 then i? S urfacc > -E-bulk- Hence, the disordered state is favorable at D < 3. On the other hand, long range order 
is favorable at D > 3. The case D = 3 is marginal and requires more quantitative consideration. We will sec that in 
this case long range order is absent for the arbitrarily weak disorder similar to D < 3. 



We begin with the analysis of the exactly solvable Larkin model |1[ . It is defined by the Hamiltonian 



H = ^ J (V<j)) 2 d 2 xdz - J h(x)<j)(x, z = 0)d 2 x (1) 

where <fi is the order parameter, h the random field, the z-axis is perpendicular to the disordered surface. At zero 
temperature we search for the energy minimum. To find the field configuration (jf{x, z) at zero temperature we 
calculate the variation of the Hamiltonian (Q) and make the Fourier transform with respect to the x coordinates: 

J(qf l -d 2 )<j ) *( % ,z) = h(q ll )6(z) (2) 

Substituting the solution of this equation 

4>*(Q\\,z) = 0*(<7||>O)exp(-|<?||k) (3) 



into Eq. (Q) we get the energy 



1 



H*= j ^[^W^M 2 -K%W{-%^)}- (4) 

From Eq. (Q) we find 4>* = cxp(— \q\\ \z)h q ., /( J\q\\ |). Assuming that the random field h is Gaussian and (h(q)h(p)} = 
AS(p + q) one finds the correlation function G{r±, r?) — {{4>{fi) — 4>{ r 2)) 2 )- For example, if z\ — Z2 <C (ri — r%) then 

G = A/(irJ 2 )ln{\x 1 -x 2 \/z 1 ) (5) 

If x\ — X2, zi <C Z2 then G = A/(2irJ 2 ) \vl{ziIz\). Since the correlation function is unlimited, there is no long range 
order. We will see that the same behavior is present in more complicated systems too, e.g. in the 3D XY model with 
surface disorder. 

As the first nontrivial example we consider a domain wall in a media with disordered surface Q . The domain wall 
is rough in the presence of bulk disorder || . We show that surface disorder also makes it rough. The shape of the 
domain wall is described by the displacement of the wall as a function of D — 1 coordinates y = </>(x, z) where x and 
z are the coordinates in a plane perpendicular to the surface, z being the direction perpendicular to the surface. The 
Hamiltonian differs from ([!]) only by the dependence of the random contribution on (j>: in the case of the random-bond 
disorder || it is a Gaussian (^-correlated random variable ^(x, (/)), (V(xi, (f)i)V(x2, ^2)) ~ 5{x\ — X2)$((f>i — We 
will demonstrate that the (4>4>) correlation function exhibits a power-law distance dependence: 



Mn)-^)) 2 )~h-^| 2( . (6) 

Introducing the field (j>* {q\ i,z) as above and repeating the derivation of Eq. (|J) we get 

f d°- 2 qu hi||</>*(gi|,0)| 2 f n , 
H * = / ,0 w5— 1" + / d D ~ 2 xV(x,cb*(x,0))- (7) 



(2ir) D - 2 2 

After replica averaging over disorder we get the effective replica Hamiltonian 

r d D - 2 % v | g ||ii^(g||,o)i 2 r , D - 3r sr WX*) - tiW) M 

Hr= J2^2^ 2f d X ^ W 2 ' (8) 

a ab 

where a and b are replica indices, R is the function describing disorder. The problem can be studied with the functional 
renormalization group (RG) following the line of Ref. g. At each RG step we integrate out the momenta from the 
interval l/(as) < q\\ < 1/a, where a is the ultraviolet cutoff and s > 1, and make the rescaling qu — > qn/s, x — > sx, 
(f>* — > s^4>*. The first term of Eq. depends on qu nonanalytically and hence does not renormalize Q. Thus, the 
temperature obeys the RG equation dT/dhiL = (3 — D — 2£)T and we find a zero-temperature fixed point. Power 
counting shows that the whole function R((f>) is a relevant operator near D = 4 in this fixed point. The RG equation 
for R can be derived exactly in the same way as in Ref. Q and to the first order in e = 4 — D has the same form as 
in Ref. § 

dR{(j))/d\nL = (e - 4QR(<f>) + C<I>R'{<P) + R"{<t>f/^ - R"{4>)R"{4> = 0) (9) 

where the factor S2/(2n) 2 = l/(27r) is absorbed into R. The solution of this equation can be read out of Ref. Q and 
gives us ( = 0.2083e in the correlation function (||). Using the result (||) for the displacement field one gets the 
same behavior (((j>(xi, z\) — (j>(x2, Z2)) 2 ) ~ \x\ ~ X2 ] \' C ' in the bulk as at the surface provided that \x\ — x?\ > z\ ~ Z2- 
Thus, the domain wall is rough in 3 dimensions in the presence of surface disorder. 

Next we consider the XY model with surface randomness. It is described by the same Hamiltonian (^) as the 
domain wall problem. However, due to the periodicity of the XY model the function R((f>) has period 2tt. One can 
apply this model for the description of vortices in a superconductor with disordered surface || . The field 4> describes 
displacements of the vortices from the regular positions in the Abrikosov lattice. Random potential is associated with 
impurities which pin vortices at the surface. 

Due to the periodicity of the Hamiltonian there is no renormalization of the order parameter, i.e. £ = 0. The RG 
equations have exactly the same structure as in the domain wall problem in dimension D + 1. Since the dimension 4 
is the critical dimension in the domain wall problem, the random XY model has the critical dimension D = 3. Thus, 
we expect the logarithmic situation in 3 dimensions: the theory is asimptotically free and R(4>) can be represented 
as R((f),L) — R* ((/>)/ In L. The exact RG equation for R* has the same structure as the one-loop RG equation (0) at 
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e = l and ( — 0. The analytical solution of this equation is known ||[| and gives the following exact result for the 
correlation function: 



(((f)(x 1 ,z 1 ) - (j}(x 2 , z 2 )) 2 ) = y lnln|aii - x 2 \ 

at |xi — x 2 \ > Z\, z 2 . The possibility to obtain an exact result makes the 3D XY model with surface disorder a good 
candidate for a quantitative test of the functional RG. Since the correlation function is unlimited at large r\ — r 2 , 
long range order is absent. However, the correlation function depends on the distance very slow. Such a situation can 
be described as super-quasi-long range order |J. 

Our treatment of the XY model in terms of the elastic Hamiltonian (||) is justified in the absence of topological 
defects. Similar to Ref. j^] we do not expect that the defects are relevant at weak disorder. At strong disorder, 
dislocations also do not proliferate deep in the bulk. Indeed, let us consider a dislocation of size L. The distance 
R of any of its points from the surface is less than the dislocation length L since otherwise the dislocation would 
not interact with the surface and should be energetically unfavorable. We estimate the energy gain associated with 
the disorder as -^disorder ~ y/AyLR, where A is the average square of the random field, and the bulk energy loss 
as £Wik ~ JL In R. One sees that the bulk loss is greater than the surface gain for large L even at strong disorder 
\fK ^> J. Thus, proliferation of large dislocation loops is unfavorable and the system renders topological order. 

In two dimensions the XY model has only quasi-long range order at low temperatures even in the absence of 
disorder. Edge disorder modifies correlation functions but does not destroy quasi-long range order as we show below. 

A recent experimental realization of a related two-dimensional system with edge randomness is a superconducting 
film with columnar defects in a part of it S. Vortices pinned by defects create a random potential near the boundary 
between the pure and impure parts of the film. This potential affects the vortex lattice in the pure part. An anisotropic 
film where vortices are free to move only in one direction can be described by the XY model with edge disorder [||. 

One cannot use the 3 — e expansion in this problem since no zero temperature fixed point exists in 2 dimensions: the 
scaling dimension of the temperature At = D — 2 becomes zero in two dimensions. Physically this result is related 
to the fact that thermal fluctuations have qualitatively the same effect on the ordering in the system as impurities: 
both destroy long-range order. Hence, we have to develop an RG procedure directly in two dimensions at nonzero 
temperature. 

Similar to the derivation of Eq. (@) we obtain the following one-dimensional replica Hamiltonian for the field <f>* at 
the edge (in the absence of vortices) 

Hr f^ m\M\ / „ EE <*°»w - a (10) 

^ ' a ab k 

where k = nk m - lu , k mm being the smallest value of k allowed by the symmetry, n an integer. Note that in this 
expression <jj* is not a zero-temperature field configuration in contrast to the previous problems. The field <f> in the 
bulk of our 2D system includes two contributions: one is related to <jf via Eq. (||) and the other is a free thermally 
fluctuating field. 

We first neglect vortices and then check how they change the behavior of the system. At T > T c = nJ/k^^ all the 
nonlinear terms in Hr are irrelevant. At high temperatures the correlation function is the same as in the absence of 
disorder and has the form ((<f>*(xi) — 4>*(x 2 )) 2 ) — 2T/(ttJ) In \xi — x 2 \. Below T c the cosine term with the minimal k 
becomes relevant. 

At fcmin = 1 we get the following one-loop RG equation 

^Ax/r 2 , T Ax 4A? 
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where a is the ultra-violet spacial cutoff. The fixed point solution of Eq. (11) is A* = (1 — T/(irJ))ir 2 J 2 /(4a). This 
result is valid, if the dimensionless combination A\a/J 2 is small. 

To obtain the correlation function we add an infinitesimal contribution A|(^^| 2 /T to the replica Hamiltonian (JlT 
The correlation function (|</>J,| 2 ) = —TdZ/dX, where Z is the replica partition function. To find the contribution to Z 
proportional to A, we calculate ^-independent corrections to Hr ([!(]) generated at each RG step. No A-dependence 
is present in the first order in Ai. Hence, we have to go to the second order in Ai. This gives us the following result: 

((<f>*(x,z = 0) - 4>* (y, z = 0)) 2 ) = — T In ^ + £ f 1 - I- Y In ^ (12) 

ttJ a 8 \ ttJ J a 
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where the first term represents the effect of thermal fluctuations and the second one represents the effect of disorder. 
At z > the correlation function has the same structure (|fj), if \x — y\ is large enough. At fc m ; n > f the result can 
be obtained from Eq. ( |H| ) by the transformation (f>* — > k(f>*, T — > k 2 T. 

The correlation function ( |l2] ) is obtained in the vortex-free model. In the absence of disorder, the Berezinsky- 
Kosterlitz-Thouless transition occurs at Tbkt = This temperature is lower than T c = itJ/k 2 ^^ in the random- 

field case fc m in = 1. Thus, vortices lead to a breakdown of our result in the presence of the random-field disorder. 
What is the effect of vortices at fc m i n > I ? To answer this question we compare the bulk elastic energy of a vortex 
in the system of size L and the energy of its interaction with the surface. The elastic energy -Ebuik ~ JhxL. An 
estimation of the disorder energy must take into account the renormalization of the disorder strength at large scales. 
The fact that we obtain a fixed point A* for A means that the effective disorder strength at scale I is of order A* a/ 1 
(rescaling at the RG steps!). The average square of the disorder energy scales hence as E 2 is ~ J dlA* /I ~ A* lni. 
There are ~ L 2 possible positions of the vortex. Assuming a Gaussian distribution for £dis we get the following 
probability of such a disorder realization that the creation of a vortex is favorable: p ~ L 2 exp(— J 2 In 2 L/ (A* In L)). 
One sees that p <C 1 , if J 2 ^S> A*. This justifies our vortex-free approximation and the result (|l2|) at the vicinity of 
the critical temperature T c as the effective disorder strength A* is small in that region. 

What happens if both surface and bulk disorder is present? This question can be easily answered in the framework 
of the Larkin model 0]. Let the average square of the bulk random field be Ab u ik- The correlation function in the 
presence of bulk disorder only jjj is ((</>(xi) — 4>{x2)) 2 ) ~ Abuikl^i — x 2 \/J 2 - This should be compared with the 
correlation function in the presence of surface disorder only (|^) . One can see that the surface random field is the main 
source of order parameter fluctuations at scales r = \x\ — x 2 \ < R c ~ A/At, u ik hi(A/(Ab u ika))- At scales r > R c bulk 
disorder dominates. 

In conclusion, we demonstrate that even arbitrarily weak surface disorder destroys long range order in systems 
of continuous symmetry at any temperature. Topological order is not destroyed and quasi-long or super-quasi-long 
range order emerges. There are two quasi-long range ordered phases in the random two-dimensional XY model. The 
results of the Letter can be relevant for the characterization of disordered surfaces because the correlation functions 
of the system in contact with the surface contain information about the surface disorder. 
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